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Periodic shedding of vortex dipoles from a moving penetrable obstacle in a 

Bose-Einstein condensate 
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We investigate vortex shedding from a moving penetrable obstacle in a highly oblate Bose-Einstein 
condensate. The penetrable obstacle is formed by a repulsive Gaussian laser beam that has the 
potential barrier height lower than the chemical potential of the condensate. The moving obstacle 
periodically generates vortex dipoles and the vortex shedding frequency fv linearly increases with 
the obstacle velocity v as fv = a(v — Vc), where Vc is a critical velocity. Based on periodic shedding 
behavior, we demonstrate deterministic generation of a single vortex dipole by applying a short 
linear sweep of a laser beam. This method will allow further controlled vortex experiments such as 
dipole-dipole collisions. 
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I. INTRODUCTION 

Quantized vortices are the characteristic excitations 
of a superfluid and they play an important role in su¬ 
perfluid dynamics [T]. Since the first realization of 
atomic Bose-Einstein condensates (BECs), vortex dy¬ 
namics has been one of the major focuses of BEC re¬ 
search [2]. In recent experiments, many efforts were 
made to study the dynamics of vortex dipoles. A 
vortex dipole is a pair of vortices of opposite circula¬ 
tions and its creation and annihilation are at the heart 
of many two-dimensional (2D) superfluid phenomena 
such as Berezinskii-Kosterlitz-Thouless superfluids 
phase transition dynamics um, and superfluid turbu¬ 
lence 13 E]. Quantized vortex dipoles were created in 
BECs by moving optical obstacles EHII] and via quench¬ 
ing dynamics nails], and also observed in quasi-2D de¬ 
generate Bose gases m- Orbital motions of a single vor¬ 
tex dipole in a trapped BEC were investigated [iQiiiiiig 
and annihilation of vortex dipoles was indirectly probed 
in relaxation of superfluid turbulence in highly oblate 
BECs [H]. 

In this paper, we present an experimental study of 
vortex shedding from a moving repulsive Gaussian laser 
beam in a highly oblate BEC. The primary result of this 
work is that when the optical obstacle is penetrable (i.e., 
the potential barrier height Vq of the laser beam is smaller 
than the chemical potential /i of the condensate), vortex 
dipoles are periodically shed from the moving obstacle. 
Furthermore, we find that the vortex shedding frequency 
linearly increases with the obstacle velocity as it exceeds 
a critical velocity. The periodic vortex shedding behavior 
is so pronounced that we demonstrate that a single vor¬ 
tex dipole can be deterministically generated by applying 
a short linear sweep of a penetrable laser beam. 

A fluid flowing past an obstacle is a textbook situa- 
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tion in fluid dynamics. In the case of a superfluid, it is 
well known that when the flow velocity exceeds a criti¬ 
cal velocity, the superfluid becomes dissipative via vor¬ 
tex generation nanHi. Early theoretical works for 2D 
geometry anticipated periodic shedding of vortex dipoles 
from a moving circular obstacle [T7H1I], where the pe¬ 
riodic shedding was depicted as a cycling process where 
the work done by a drag force from the moving obstacle 
accumulates to generate a quantized vortex dipole of a 
finite energy. However, in previous BEC experiments 0- 
ini, periodic vortex shedding has never been observed, 
although the existence of a critical velocity was clearly 
demonstrated. 

In generating a vortex dipole, a penetrable obstacle 
with Vq/ jii < 1 has a clear advantage over an impenetra¬ 
ble obstacle. Because there is no density-depleted region 
in the condensate, it is ensured that a vortex and an an¬ 
tivortex would be created at the same time into the con¬ 
densate. On the other hand, in the impenetrable case, 
there is a density-depleted region occupied by the obsta¬ 
cle such that a vortex can be individually emitted into 
the condensate, while the other vortex that is created 
as a partner of the emitted vortex still remains in the 
density-depleted region [22l |23] . The individual vortex 
emission might be induced by small fluctuations in su¬ 
perfluid flow, braking the mirror symmetry of the system 
with respect to the obstacle moving direction. Recent nu¬ 
merical studies on vortex shedding from an impenetrable 
obstacle showed that vortex shedding is mostly irregular 
and regular shedding of vortex dipoles occurs only in a 
stringent condition of the obstacle [22j [24] , which is prac¬ 
tically difficult to satisfy in current experiments. In this 
work, we observe that a moving penetrable obstacle sheds 
vortex dipoles in a periodic and stable manner. 

This paper is organized as follows. Sec. II describes 
our setup and procedure of the vortex shedding exper¬ 
iment. Sec. Ill presents the measurement results for 
various moving conditions of the obstacle and discusses 
the difference between the penetrable and impenetrable 
cases. Sec. IV demonstrates an experimental method for 
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FIG. 1: (Color online) (a) Schematic of the experiment. A 
penetrable obstacle is formed by a repulsive Gaussian laser 
beam and it moves with velocity v in the center region of a 
highly oblate Bose-Einstein condensate (BEG). The barrier 
height Vq of the optical potential is lower than the chemical 
potential fi of the BEG. (b) Time-of-flight absorption image 
of a BEG. A vortex dipole is generated by the moving obstacle 
and it is identified with its pair of density-depleted cores. 


deterministically generating a single vortex dipole and 
examines its stability in terms of of the linear velocity of 
the generated vortex dipole. Finally, Sec. V provides a 
summary and outlook. 


II. EXPERIMENT 

The experimental setup is similar to that in our pre¬ 
vious work m- A highly oblate BEC of ^^Na atoms is 
prepared in a hybrid trap composed of a pancake-shaped 
optical trap and a magnetic quadruple trap. The trap¬ 
ping frequencies are oor^z = 27r x (9.0,400) Hz and the 
atom number in the condensate is Nq = 3.3(2) x 10^, 
giving the chemical potential /i ~ x 1.0 kHz, where h is 
the Planck constant. The condensate fraction of the sam¬ 
ple is over 80%. At the trap center, the healing length 
is ^ = fi/y/2m/i 0.45 jim and the speed of sound is 
Cs = \//V^ ~ 4.3 mm/s, where h = h/27r and m is the 
atomic mass. The T homas-F ermi radius of the trapped 
condensate is R = 110 /im. 

An optical obstacle is formed by a repulsive Gaussian 
laser beam whose 1/e^ width is cr = 9.1(12) jim ^ 20^. 
We initially locate the laser beam 11.5 /im away from the 
condensate center and horizontally translate it through 
the center region of the condensate using a piezo-driven 
mirror [Fig. 1(a)]. The moving velocity v of the laser 
beam is kept constant during the translation. After the 
laser beam sweep, we linearly ramp down the beam power 
within 20 ms and then take an absorption image of the 
condensate to detect vortices [Fig. 1(b)]. When we re¬ 
lease the trapping potential, we turn off the magnetic 
trap 13 ms earlier than the optical trap and place 24 ms 
free expansion before taking the image. 


III. PERIODIC VORTEX SHEDDING 
A. Penetrable obstacle 

We first investigate the case of a penetrable obstacle 
with Vo/jii ~ 0.74 by measuring the number of vortex 
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FIG. 2: (Golor online) Periodic shedding of vortex dipoles 
from a moving penetrable obstacle of Vb//^ ~ 0.74. (a) Num¬ 
ber of generated vortex dipoles, Nd as a function of the sweep¬ 
ing time t for various moving velocities v of the obstacle. Each 
data point consists of 15 realizations of the same experiment 
and the error bar indicates the standard deviation of the mea¬ 
surements. The activation time Tv for generating the hrst 
vortex dipole was determined from a sigmoidal function ht 
(red solid lines, see the text for detail) to the data up to the 
Nd=l plateau region. The shaded area indicates the region 
of < t < 2tv . (b) Vortex shedding frequency R = 1/tv SiS 
a function of v. The solid line is a linear ht to the data. 


dipoles generated from the obstacle as a function of the 
sweeping time t for various moving velocities v. The mea¬ 
surements results are displayed in Fig. 2(a). We see that 
the vortex dipole number shows stepwise increasing 
behavior as a function of t. There is a certain amount 
of an inert time Ty before generating the hrst vortex 
dipole and another sweeping time almost equal to Ty is re¬ 
quired to generate the second vortex dipole. This regular 
stepping-up feature clearly indicates that the inert time 
is not due to the transient response of the condensate 
to the perturbations of the laser beam but an intrinsic 
activation time for generating a vortex dipole from the 
moving obstacle. The periodic vortex generation might 
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FIG. 3: (Color online) (a) Vortex dipole number Nd as a 
function of v for fixed sweep length L = vt. Vo/m ~ 0.74. The 
blue solid lines are the theoretical predictions obtained from 
a model assuming perfect periodic vortex shedding with fv = 
a{v — I’c), a = 0.25 /im“^, and Vc = 1.0 mm/s [Fig. 2(b)]. (b) 
Critical velocity Vc measured for fixed L. Vc was determined 
from the data in (a), using the data analysis method described 
in Ref. m The solid line is a linear fit to the data. 


be understood as a cycling process where the moving ob¬ 
stacle exerts a drag force during the activation time and 
the accumulated energy is released as a vortex dipole at 
a certain threshold condition m- 

It was theoretically anticipated that the vortex shed¬ 
ding frequency fy linearly increases with the moving ve¬ 
locity V of the obstacle as fy = a{v — Vy), where Vy is 
a critical velocity [m [HHH113 [IS- Interpreting the 
activation time Ty as the inverse of fy, we find that 
our experimental data support the theoretical predic¬ 
tion [Fig. 2(b)]. Ty was determined by fitting a sigmoidal 
function = [1 + to the data up to the 

Nd=l plateau region. Here 7 represents the jittering of 
the vortex emission event and we used the value of 1.67 
as the measurement uncertainty of Ty. Fig. 2(b) shows 


fy = l/Ty as a function of v and the result is well de¬ 
scribed with the functional form of fy = a{v — Vy), giving 
Vy = 0.99(4) mm/s and a = 0.25(2) /im“^. 

We present another set of measurements in Fig. 3(a), 
where is measured as a function of v for fixed sweep 
lengths L = vt of the laser beam. This is a typical situ¬ 
ation addressed in previous experiments to measure the 
critical velocity Vy for vortex shedding unimiis]. Our 
motivation is to examine the effect of the finite activa¬ 
tion time Ty on the determination of Vy. In order to 
generate vortices, the sweeping time t = Ljv should be 
longer than Ty, which requires v < aL{v — Vy) from the 
relation fy = a{v — Vy) = l/ry. Therefore, in the fixed 
L setting, the critical velocity would be measured with 
a systematic error of 5vy = VyjiaL — 1). Indeed, we 
observe that the measured critical velocity Vy increases 
with decreasing L [Fig. 3(b)]. Vy was determined us¬ 
ing the data analysis method described in our previous 
work m- The inhomogeneous density effect is negligible 
because the variation of the local speed of sound Cg is 
less than 5% even for the longest L = 46 /im. The linear 
extrapolation of the experiment data to l/L = 0 gives a 
velocity of 0.97(3) mm/s [Fig. 3(b)], which is remarkable 
in that this value is consistent with the critical velocity 
Vy obtained from the vortex shedding frequency data in 
Fig. 2(b). 

The blue solid lines in Fig. 3(a) indicate the model 
predictions assuming perfect periodic generation with 
fy = a{v — Vy), a = 0.25 /im“^ and Vy = 1 mm/s. Most 
of the qualitative features of the experimental data are 
explained with the model, including the reduction of the 
Nd=l plateau region with increasing L. This observation 
corroborates the periodic shedding of vortex dipoles from 
the moving penetrable obstacle. 


B. Impenetrable obstacle 

We perform the same measurements with an impene¬ 
trable laser beam of Vb//^ ~ 1-5 (Fig. 4). Interestingly, 
the activation time is still necessary for generating the 
first vortex dipole and reduces with higher moving veloc¬ 
ity as observed with the penetrable obstacle. This seems 
to be consistent with our description that the activation 
time is required to accumulate a certain amount of en¬ 
ergy with the drag force to form a vortex dipole. How¬ 
ever, no periodic vortex shedding is observed with the 
impenetrable obstacle and the number of vortex dipoles 
grows monotonically once the first vortex dipole is emit¬ 
ted. This means that the subsequent vortex generations 
are significantly affected by the presence of the first vor¬ 
tex dipole [23l|27]. It is noted that for v = 1.90 mm/s in 
Fig. 4(a), the increasing rate of is almost four times 
higher than the inverse of the activation time, which im¬ 
plies that vortex shedding can be accelerated in the pres¬ 
ence of other vortices near the obstacle. 

The distinctive difference of the impenetrable case 
from the penetrable case is that there is a density- 
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FIG. 4: Vortex shedding from a moving impenetrable obstacle 
of Vo/m ~ 1-5. (a) Vortex dipole number Nd versus sweeping 
time t. (b) Nd as a function of the moving velocity v for sweep 
length L = 23 /im. 

depleted region occupied by the obstacle. Hence, it is 
possible for vortices to be emitted individually into the 
condensate. The individually emitted vortices would dis¬ 
turb the superfluid velocity field around the obstacle 
more strongly than a vortex dipole, leading to irregu¬ 
lar vortex shedding or possibly more intriguing shed¬ 
ding patterns such as the Benard-von Karman vortex 
street [22] [23], 

In Fig. 5, we present several vortex images containing a 
few vortex dipoles. It is clearly seen that the vortex shed¬ 
ding pattern is irregular in the impenetrable case and was 
often observed that the propagation direction of the first 
vortex dipole largely deviated from the moving direction 
of the obstacle [Fig. 5(a)], implying that the two vortices 
were emitted at different times from the obstacle. Fur¬ 
thermore, vortex clusters having larger density-depleted 
cores were sometimes observed [Fig. 5(c)] [101 EH- Oii 
the other hand, the first vortex dipole from a penetrable 
obstacle usually propagates along the moving direction of 
the obstacle. Symmetric periodic array of vortex dipoles 
are known to be unstable towards staggered configura¬ 
tions of vortex dipoles [22l |23l [28] , which were typically 
observed in the experiment when multiple vortex dipoles 
were generated [Figs. 5(d) and (e)]. Occasionally, sym¬ 
metric trains of vortex dipoles were observed [Fig. 5(f)]. 

IV. DETERMINISTIC GENERATION OF A 
SINGLE VORTEX DIPOLE 

In the experiment data for the penetrable obstacle in 
Fig. 2(a), it is evident that fluctuations of the vortex 
dipole number are strongly suppressed for Ty < t < 
2ry and even vanish for some conditions. This is an out¬ 



FIG. 5: Various patterns of vortex shedding from (a)-(c) im¬ 
penetrable and (d)-(f) penetrable obstacles. In the impen¬ 
etrable case, the first vortex dipole was often observed to 
propagate in a direction largely deviated from the -\-x mov¬ 
ing direction of the obstacle. The dashed lines indicate the 
moving line of the obstacle, (c) Vortex clusters having larger 
density-depleted regions were also observed. In the penetra¬ 
ble case, staggered configurations of multiple vortex dipoles 
were typically observed, (f) Occasionally, symmetric trains of 
vortex dipoles also appeared. 


standing feature of the periodic vortex shedding, suggest¬ 
ing that deterministic generation of a single vortex dipole 
is possible when applying a short sweep of the penetra¬ 
ble laser beam. In this section, we characterize this short 
sweep method by examining the stability of the propa¬ 
gation velocity of the generated vortex dipole. 

In Figs. 6(a) and (b), we show two exemplary images 
for V = 1.52 mm/s with t = 11.5 ms and 7.5 ms, re¬ 
spectively. In the first image, a vortex dipole is located 
near the center region and clearly identified with its two 
distinctive density-depleted cores. On the other hand, 
in the second image that is taken with a short sweep¬ 
ing time t ^ Ty^ a small crescent-shaped density dimple 
is observed with low visibility in the boundary region of 
the condensate. We interpret this density dimple as an 
incipient vortex dipole with a short intervortex distance 
d, which cannot be fully developed because of the short 
sweeping time. The incipient vortex dipole moves fur¬ 
ther away from the sweeping region because of higher 
propagation speed Vd = h/md for smaller d and its low 
visibility might be attributed to the blurring effect during 
the expansion before imaging. 

Compiling all the images containing a single vortex 
dipole for v = 1.52 mm/s, in Fig. 6(c) we display the 
distribution of the vortex dipoles in the plane of radial 
position and visibility. The visibility was determined as 
1 — n^/n, where rim is the density minimum value in the 
core region and n is the average density in the surround¬ 
ing area of the vortex dipole. The distribution shows 
a strong correlation between the visibility of the vortex 
dipole and its radial position, i.e., propagation speed, re¬ 
vealing the temporal development of the vortex dipole. 
From a rough estimation, taking into account the ex¬ 
pansion time in the imaging sequence, we find that the 
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FIG. 6: (Color online) Formation of a vortex dipole. Images 
of BECs containing a vortex dipole, where Vq/m ~ 0.74 and 
the sweeping times were (a) t = 11.5 ms and (b) 7.5 ms at 
V = 1.52 mm/s, respectively. For long t, the vortex dipole 
is clearly identified with its two vortex cores near the center 
region. For short t, it appears as a small density dimple with 
low visibility in the boundary region. The enlarged images 
of the vortex dipoles are presented on the right side of each 
image, (c) Vortex distribution in the plane of radial position 
and visibility, obtained from all the images containing a single 
vortex dipole in the experiment for v — 1.52 mm/s [Fig. 2(a)]. 


propagation speed of the most incipient vortex dipole is 
comparable to the moving velocity of the obstacle, which 
is expected from the fact that the vortex dipole is shed 
from the moving obstacle. The condition oi Vd = v gives 
the intervortex distance d = h/mv ^ 4^. 

When the sweeping time is not sufficiently long, al¬ 
though long enough to ensure generation of a vortex 
dipole, the speed of the generated vortex dipole is not 
reproducible. For example, when t = 8.5 or 10 ms, 
the vortex generation probability is almost unity but 
the vortex dipole velocity widely scatters [open circles in 
Fig. 6 (c)]. This is partially due to the uncertainty of the 
vortex emission time and also probably due to the finite 
dwelling time of the vortex dipole in the proximity of the 
laser beam, which might be estimated to be At = a/2v. 
Therefore, for deterministic generation of a single vortex 
dipole, the sweeping time needs to be optimally chosen 
to have better stability in the vortex dipole velocity Vd- 
We examined the stability of Vd by putting an additional 
100 ms hold time before releasing the trapping potential 
and measuring the change in the vortex position. In our 
experiment, the propagation speed of the vortex dipole 
was stable within 10 % and its propagation direction was 
parallel within 15 degrees to the sweeping direction of 
the laser beam. 

It needs to be mentioned that we cannot exclude the 
possibility that the crescent-shaped density dimple is a 
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FIG. 7: Dependence of the intervortex distance d on various 
sweeping parameters of the penetrable obstacle, (a) d versus 
V, where a = 9.1(12) pm and Vb/p ~ 0.74. (b) d versus 

Vb/p, where a = 10.7(13) pm and v = 1.25 mm/s. (d) d 
versus cr, where Vb/p ~ 0.74. In (d), v was controlled to 
range from 1.08 mm/s to 1.34 mm/s due to the change of the 
critical velocity m- In each measurement, the sweeping time 
t was chosen to have d stable (see the text for detail), d was 
measured in the image and divided by the scaling factor of 
2.1 for the expansion in the imaging sequence. 


rarefaction pulse or a gray soliton, which can also prop¬ 
agate with preserving its density profile with a nonzero 
density minimum [25j |29]-[3T] . When the energy accumu¬ 
lated by the moving obstacle is not sufficient to generate 
a vortex dipole, it can possibly transform into lower en¬ 
ergy excitation. The sudden stop of the moving obstacle 
might play a role in forming and stabilizing the excita¬ 
tion. 

Finally, with a motivation to find a way to control the 
velocity of the vortex dipole, we investigate the depen¬ 
dence of the intervortex distance d on various parameters 
of the moving obstacle. Fig. 7 displays the measurement 
results of d as functions of the moving velocity the 
barrier height Vb/p, and the beam waist a of the laser 
beam. We see that d slightly decreases with increasing 
suggesting one way to control Vd but in a limited range. 
Intriguingly, d is almost insensitive to Vb/p and a. Con¬ 
sidering the fact that the critical velocity Vc is signifi¬ 
cantly affected by the potential parameters of the obsta¬ 
cle m, the observed weak dependence of d on them is 
quite unexpected. In Ref. m, vortex shedding from a 
moving attractive obstacle of Vb/p < 0 was numerically 
studied, showing that d increases with higher which is 
opposite to our observation. Further theoretical studies 
are warranted to understand what determines the size 
of the vortex dipole when it is generated from a moving 
penetrable obstacle. 


V. SUMMARY AND OUTLOOK 

We investigated vortex shedding dynamics with a mov¬ 
ing penetrable obstacle in a highly oblate BEC. We ob¬ 
served periodic shedding of vortex dipoles and found the 
linear dependence of the vortex shedding frequency on 
the moving velocity of the obstacle above the critical ve¬ 
locity. In addition, we demonstrated deterministic gener¬ 
ation of a single vortex dipole by applying a short linear 















6 


sweep of the penetrable laser beam. This vortex gener¬ 
ation method is expected to provide new opportunities 
for further controlled experiments on vortex dynamics. 
For example, as suggested in Ref. [30], if two laser beams 
are employed to generate two vortex dipoles separately 
at different positions in a condensate, it would be pos¬ 
sible to investigate collision dynamics of vortex dipoles 
in a controlled manner. Dipole-dipole collisions are par¬ 
ticularly interesting in that vortex pair annihilation may 
occur during the collision. Vortex pair annihilation is 
one of the key issues in the study of 2D quantum tur¬ 
bulence [321440] . Many aspects of the pair annihilation 


process mnn], such as the roles of vortex-phonon inter¬ 
actions Uni [44] , and thermal dissipation [45] , have never 
been directly addressed in experiments. 
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